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Results are presented for the nonequilibrium infinite-U Anderson model using a large-N approach, 

where N is the degeneracy of the impurity level, and where nonequilibrium is established by couphng 

the level to two leads at two different chemical potentials so that there is current flow. A slave-boson 

representation combined with Keldysh functional integral methods is employed. Expressions for the 

static spin susceptibility xs and the conductance G are presented to O (■^) and for an applied 

'q«^' ■ voltage difference V less than the Kondo temperature. The correlation function for the slave-boson 

' is found to be significantly modified from its equihbrium form in that it acquires a rapid decay in 

, time with a rate that equals the current induced decoherence rate. Physical observables are found 

. to have a rather complex dependence on the coupling strength to the two leads which can lead 

' to asymmetric behavior XsiV) 7^ Xs{—V), G{V) 7^ G{—V) both in the mixed valence and in the 

M ' Kondo regime. 

P ■ 

^ : I. INTRODUCTION 

The theoretical problem of strong correlations coupled with nonequilibrium has become an active area of research in 
' recent years, in part due to the enormous success in realizing experimental systems which can be driven out of equilib- 
rium in a controlled manner. Some examples of these are current carrying quantum dots and single molecule devices^, 
strongly driven ferromagnetic systems^,, cold atoms trapped in optical lattices with rapidly tunable parameters^. One 
of the theoretical challenges in the study of out of equilibrium strongly correlated systems is that, unlike systems in 
• equilibrium which are characterized by some underlying principles such as the energy minimization principle, no basic 
^ underlying principles are known for out of equilibrium systems making it rather difficult to develop general theoretical 
I , techniques to study them. 
' O ' Perhaps the most actively studied out of equilibrium systems are nonequilibrium quantum impurity models which 
C , are systems characterized by a few local degrees of freedom coupled to one or more reservoirs (as in a quantum dot 
or a molecular conductor), and where nonequilibrium is achieved by maintaining the reservoirs at different chemical 
potentials and/or by subjecting the system to time-dependent fields. For strong local interactions the ground state 
of quantum impurity models show many-body resonances such as the Kondo or polaronic resonance. The effect of 
current flow on these resonances has been studied using a variety of methods such as renormalized perturbation theorjsi, 
, flow equation methods^, real time renormalization group on the Keldysh contour— il, and functional renormalization 
' group methods*. While these approaches are applicable when the external drive is large as compared to the Kondo 
, temperature, in the opposite limit of drive small compared to the Kondo temperature, perturbative methods based 
^—1 ' on Fermi- liquid theory have been uscd^. There have also been efforts at developing exact solutions based on the 
' construction of exact scattering states in the presence of current flowiSiii. There are also several promising numerical 
methods that are being developed such as the real-time numerical renormalization group method^^, quantum Monte 
^\ Carlo computation of real time Keldysh diagramsi^^i^, iterative summation of real time path integrals^^ and the 

imaginary time formulation of real-time nonequilibrium problems 
^ In this paper we will use large-N methods^^ to study a nonequilibrium quantum impurity model. In particular, 

• ^ , we will study the Anderson model when the on-site Coulomb interaction U = 00, and in addition the system has 
■ been driven out of equilibrium due to current flow. N here will represent the degeneracy of the impurity level. The 
J-j _ physical systems this corresponds to are quantum dots or molecular devices where the level active in transport is 
■ - - • characterized by a total angular momentum J = L + S which is large, and hence has a large degeneracy N = 2J + 1. 

This could arise due to the particular form of the confining potential in the quantum dot, or by the use of a molecule 
where conduction occurs via a metal ion with a partially empty outermost d or / orbital. Note that the infinite-U 
Anderson model under out of equilibrium conditions has so far been studied using the non-crossing approximation 
(NCA)l* and slave boson mean-field methods — . In this paper we will also employ the slave-boson representation 
which is a convenient way to project out all states except the empty and singly occupied state of the dot^°. However, 
we will go beyond mean-field by including the effect of fluctuations to O (jj)- Our theoretical approach is closest to 
that of Read et al.-^'^^ , but carried out for a nonequilibrium system using Keldysh functional integral methods. 

A few words on the regime of validity of the results presented in this paper. The U — 00 limit of the Anderson 
model is the so called mixed-valence regime where the system is characterized by both local charge as well as spin 
fluctuations. The Kondo regime may be accessed by making the bare level energy large and negative in which case the 
charge fluctuations arc frozen out and only the spin-fluctuations exist. This limit can be taken in a straightforward 
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way in all physical observables. Thus we will present results for the nonequilibrium static spin susceptibility and the 
conductance in both the mixed valence as well as in the Kondo regime. 

This paper is organized as follows. In section we present the model, and write it as a Keldysh path-integral 
suitable for studying nonequilibrium systems. In section IIIII we briefly present the main results of the paper before 
turning to the full calculation. In section IIVI we study the Keldysh path integral in the limit of iV — + oo when a 
mean- field or saddle-point approximation becomes exact. In this limit the voltage dependence of the local charge 
density, static susceptibility and the conductance are derived. Following this, the rest of the paper is devoted to 
the study of the effect of fluctuations to 0{\/N). As found by Read et a^, the 1/N corrections are in general 
associated with infra-red divergences whose origin is the zero-mode of the slave-boson representation. While the 
infrared divergences are logarithmic in equilibrium, we find that out of equilibrium the divergences become more 
severe with a pole structure. However, just as in equilibrium, in the computation of all physical observables these 
infrared divergences are found to cancel so that the final expressions are well defined. 

The 0{1/N) computation is organized as follows. In section |V] the mean- field saddle point expressions for the 
level position and the level broadening are corrected to 0{1/N). In Section IVTl the local impurity charge density is 
computed. In Section [VIII the bosonic correlation fimction is evaluated. While in equilibrium the bosonic correlation 
function has a power-law decay in time with an exponent consistent with X-ray edge physics^^, for the current carrying 
case we find that the long time behavior has both a power-law as well as a rapid exponential decay in time, the latter 
arising due to current induced decoherence. The bosonic correlation function appears in the computation of various 
physical observables. We present results for the static susceptibility in section IVIIIl while expressions for the impurity 
spectral density and conductance are presented in IIXI Many of the details of the computation are relegated to the 
appendices. Finally we conclude in section 1X1 



II. MODEL 



We use the slave-boson representation— of the infinite-U Anderson model which is a convenient way to project out 
all except the empty and singly occupied states of the impurity level. The Hamiltonian in this representation is, 



(1) 



k,m.a — L^R 



where m = — J ... J represents the spin-projection of the local level, N=2J-|-1 is the degeneracy of the level, Ckma 
represent the lead electrons, and we have generalized to the case where there are two leads (labeled by a = L, R) 
which will be maintained at two different chemical potentials ^il.r to capture the nonequilibrium current carrying 
case. is the hybridization to the two the leads. The above Hamiltonian is accompanied by the constraint 



(2) 



to ensure that the system remains within the restricted Hilbert space of an empty or singly occupied local level. 

We write the Keldysh path integral^^^ for Eq. [1] and impose the constraint in Eq. [2] by introducing two Lagrange 
multipliers A± 

(3) 



Zk — / 2? [/m±, /m±, A±, 6±, 6^, Cm, Cm] exp (iTr[SK]) 

where the Tr symbol in Eq. [3] represents a trace over time indices, and 
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In the above g^^ is the inverse Green's function for the leads and is a 2 x 2 matrix in Keldysh space. It is convenient 
to integrate out the lead electrons to obtain, 
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Performing a rotation to retarded (R), advanced (A), Keldysh (K) space^'^, and defining the quantum fields as 
Oq = (0_ — 0+)/2 and the classical field as Od = (0_ + 0+)/2, we get 



SK—^'^{frn,q fm.cl) f — XcITq — XgTx — {bdTQ + bgTx) 'Sc {b*ciTo + b*Tx) 
m 

where the Ec are the self-energies due to coupling to leads, 

^'^ ~ ' E;* 



fm,cl 



(6) 



(7) 



with J:i="'^'^{t,t') = i Y.k,a=L,RVa9i""'^'^ik;t,t'). Thus the self-energies due to coupling to leads is O (^). We 
will make the assumption of constant density of states in the leads which gives 



a=L,R 

Sf = -2*r ^ £fi (1-2/(0; -Ma)) 

a=L,R 



(8) 

(9) 
(10) 



The aim will be to use the action in Eq. [5]to evaluate physical observables perturbatively in Before turning 

to the full computation, we present the main results in the next section. 



III. BRIEF DISCUSSION OF RESULTS 



Let us suppose that the chemical potential of the left lead is /il = V/'^ while that of the right lead is fin = —V/2. 
As specified in Eq[8l let rL(rfl) be the self-energy due to couphng to the left (right) lead, while F = -f F/j is the 
total self-energy. In terms of the above parameters, the static susceptibility in the Kondo regime (denoted by the 
superscript = 1 to indicate the value of the charge on the impurity level) is found to have the following universal 
form, 
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(11) 



sing the mean-field 



Kondo temperature — , and the Csi are numbers specified in the text (after Eq. I146P . Thus one finds that for an 
asymmetric couphng to leads (Fl ^ Tr), XsiV) 7^ Xsi—V)- This lack of symmetry when V ^ —V arises due to the 
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fermi-level dependence of the Kondo temperature. To see this we set the couphng to one of the leads (say Tji) to 
zero. This corresponds to an equilibrium configuration where there is no current flow. For this case Eq. [11] reduces to 



xT 



--\fiL = V/2,TR^0) 



g^filJ{J + l) 



3T, 
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N \2Tk 
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(12) 



Thus the terms in Eq. [12] can be interpreted as a change in the Kondo temperature arising from a change in the 
chemical potential of the left lead by (5/ii, = V/2. The asymmetry xsiV) 7^ Xs(— ^) in the Kondo regime therefore 
arises when the level is unequally coupled to two leads, each associated with a different equilibrium Kondo temperature. 

In contrast, the terms of the type ^p^" (^) Eq[TT]are purely nonequilibrium terms that arise due to inelastic 
scattering processes in the energy window V when there is current flow, and are thus associated with current induced 
decoherence. The identification of these terms with decoherence becomes clearer below when we discuss the slave- 
boson correlation function. 

We now turn to the discussion of the conductance. Here too one finds that the fermi-level dependence of the spectral 
density can give rise to a conductance that is asymmetric under V — V^^. In particular the mean-field saddle point 
expression for the conductance in the mixed valence regime is found to be 



Gsp{V) = Gsp{V = 0) 
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(13) 



where JT-i? is the charge density on the level when — fJ-R — 0, and GgpiV = 0) ^^^f^ i.^^jf'Y ■ "^^'^ conductance 
in the Kondo regime can be accessed by taking the limit n^? — > 1 in Eq. 1131 Thus for a symmetric coupling to the 
two leads, the mean-field conductance in the Kondo regime becomes 



G--=\v- Tl = r^) = Gir\v = 0) 



1 - 



V 



(14) 

for the mixed 



The 1/iV correction to the conductance for the case of symmetric couplings to leads is given in Eq, 
valence regime and in Eq. 11661 in the Kondo regime. 

We now turn to the discussion of the bosonic correlation function DK{t,t') = ~i{{b{t),b^ {t')}) which is used to 
obtain the physical observables discussed above. At the mean-field level, b{t) (6) in the Hamiltonian (Eq.[T]), so that 
the U{1) symmetry of the Hamiltonian is broken. In equilibrium, including fiuctuations to 0{1/N), the correlation 
function become a^^'^^ 



D%{t) = ~2^{l-nF){l-^\n{tn] 



(15) 



It was argued that^^'^^ since the model in Eq.[l] cannot have any broken symmetry state, including terms to higher 
orders in 1/A^ should lead to a power-law decay in the bosonic correlation function so that the symmetry of the 
Hamiltonian is restored. Thus, 



Dm 



[tTir 



(16) 



where a 



N 



and equals the Nozieres-de Dominicis infrared exponent for the response of an electron gas 
subjected to a sudden change in potential^^. 

We find that the result for the bosonic correlation function for the current carrying case, and for long times {Vt 3> 1) 

is, 



where the coefficients cl,r,, 
terms of O {V/Tfj are. 
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are specified in Eqns. 11251 11261 The Ci are weakly voltage dependent, and neglecting 
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For evaluating quantities to 0{1/N), Eq. [T7]is sufficient. However it is interesting to consider how would change 
when higher order in 1/N terms are included. Following Eq. 1161 we expect that the bosonic correlation function will 
have the form 

- -P (~f ^^*) (20) 

where aneq = (cL+Cfj) /N. Thus to all orders in 1 /N the bosonic correlation function will be characterized with a long 
time power-law decay along with rapid exponential decay in time, the latter arising due to current induced decoherence. 
The rate of decoherence is '"p^" V, and is a energy scale that appears repeatedly in all physical observables. Note that 
Eq. [20] is also consistent with nonequilibrium X-ray edge physics i.e., the response of an out of equilibrium electron 
gas to a sudden change in potential studied recently in various contexts^^. 
We now turn to the derivation of the above results. 



IV. MEAN-FIELD SADDLE POINT TREATMENT 



In the mean field saddle point treatment, one assumes the fields hci,q,Xci,q in Eq. [6]to be constants in time. The 
action Sk is then minimized both with respect to the classical fields 5cz, \ci and the quantum fields bq, Xq. The classical 



saddle points — 0, — are automatically satisfied for bq = Xq = 0. Thus in order to satisfy the saddle point 
equations with respect to the quantum fields = 0, = it is sufficient to expand Sk to linear order in the 
quantum field. To carry these steps out, we integrate out the fermionic fields in Eq. [6]to obtain 



-iNTr In 



^nif ~ ^q'^x ~ bcl^cbqTx — b qT,j.Yi J:) d 
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where the mean-field fermionic Green's function is 



Defining 



where F plays the role of the level broadening. 
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From Eq. [2T1 the saddle point equation for Xq, — gives 
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Using Eq. [^nithe above becomes. 
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After performing the frequency integrations, we obtain 



1 = 



r N 



r TT 

Similarly, minimizing Eg. 1211 with respect to bq ^i leads to 
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Using expressions for Sc, the above leads to 



(29) 



(30) 
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which after performing the frequency integrations gives. 
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We will now proceed to solve the two saddle point equations Eq. [51] and [321 and use the solution to evaluate various 
observables. The results obtained will be exact in the limit N ^ oo. 

Solution of the saddle point equations 
Let us define 



Ac/ + Eq ^ ep 



(33) 



where ep is the effective position of the impurity level. When N — > oo, F, F ^ and iVF = const. Using this, Eq. 
may be simplified to 



F iVF 



Fi/F , Tr/T 



while Eq [321 becomes (defining ep ~ Tj^ as the position of the level in the limit N oo) 
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(36) 
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m should not to be confused with the label for the spin projection. Note that in equilibrium, ml = M-R = Oi '"^v = "m. 
In terms of these variables, Eq. [331 implies the following for the saddle point solution for bd 



"^'"^^f " F - l + mv 



whereas the impurity charge density is 



up = --^^Tr [G^f] + ^ ^ 1 ^ 



my 



2 ^ "'^J 2 
Note that in the Kondo limit, my ^ 1 so that np ^ 1. 



F 1 + my 



(38) 



(39) 
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A. Solution for Ta 



We solve Eq. [35] when -Eq > Ta- Writing 



TA = n+ STa 



(40) 



where 



= De — ~ 



(41) 



is the equilibrium solution for the impurity level, Eq. [35] becomes 

in + STa - A^Lr^/^m + sta ~ M/^r^/^ = n 



(42) 



For small voltages, |(5Ta|, |/iL,fl| << T^, a Taylor expansion leads to the following expression for the change in Ta 
due to bias, 



(43) 



B. Mean field impurity susceptibility 



We now turn to the evaluation of the voltage dependence of the impurity susceptibility. The spin-response function 
at the mean-field level is given by 



Am/ 



(44) 



Using the identity X]m=-,/...,/ 
spin-response function becomes. 



-^J{J + 1) = ^J{J + 1), the spin susceptibility which is the zero frequency 



Xs, = xU^ = 0) - ^ J(J + 1) 



9 ? 



(TA-ML)' + f2 {TA-mf + '^^ 



(45) 



For ^ oo we may drop terms of 0{T^), 



Xs, = ^-^J{J+1) 
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(1 + mv)TA 



r./r 



(46) 



Taylor expanding Eq. [46l in powers of ^T^fr- and defining 



mo 



NT 



we find the following voltage dependence of the susceptibility at saddle-point, 



(47) 
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In the Kondo limit, uiq ^ 1, or the equilibrium charge on the level np — 
susceptibility becomes 



(48) 

1. In this case the static 
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Mean-field conductance 



The current is given bj*^ 

ie 2TlTr 



I = 



r 



where 



Within mean-field, 6± are constants in time and equal to the saddle point value given in Eq. 
temperature Eg. [501 becomes. 
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(50) 
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Thus at zero 
(52) 



Let us set hl = e.V/2,^ii = —eV/2. The zero-bias conductance depends only on the equilibrium properties of the 
spectral density and is given by, 
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(53) 



The non-linearity in the conductance arises due to the frequency and voltage dependence of the spectral density 
(namely the voltage dependence of its position Ta and its width F). We find the following expression for the non- 
linear conductance. 
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(54) 



The above implies that for asymmetric coupling to the leads (F^ 7^ Fij), the conductance shows a rectification type 
behavior, i.e. Gsp{V) ^ Gsp{—V). Whereas for symmetric couplings to the leads, the conductance reduces to 



Gsp{V;TL = rR)^Gsp{V = 0) 



3 mo 



V 



The conductance in the Kondo limit can be obtained by setting tuq 3> 1 in Eq. 
we get 



G" 



V_ 

4 vn 



(55) 

Thus for symmetric couplings, 

(56) 



The main results of this section are the expressions for the static susceptibility (Eq. 251 H^ . and the conductance 
f Egns. [54l [55l [56)1 . In the rest of the paper we will study how these results are modified when fluctuations to O (■^) 
are taken into account. 



V. FLUCTUATIONS ABOUT MEAN-FIELD 



We now turn to the computation of how the saddle-point Eqns [291 and |32l get modified when fluctuations are 
included. Formally the steps involved are to write bd bgp -f bd, 6*; bsp -f 6*;, bq bq + bq, b* — > bq + b*. 
Then we integrate out all the fermionic and bosonic fields bd , 6*; ,bq,b* obtaining a resulting action that depends 
only on Sk = Sxib^p^bq, Xd, ^q)- Each of the variables x = bsp,bq, Xd, Xq are then determined by requiring that 



FIG. 1: Diagrams representing the mean-field Green's function Gmf, the self-energy due to coupling to leads Ec, and the 
bosonic propagator D. Each has a 2 x 2 Keldysh structure. 



— g^^" = 0. Of course, the bosonic and fermionic fields cannot be integrated out exactly. This is therefore done 
perturbatively in j^. Moreover, as discussed in Section HVl the saddle point equations with respect to the classical 
fields — 0, 1^ — are always satisfied if all the quantum fields \q ~bq~ 0. Thus to obtain the quantum saddle- 
points, it suffices to expand Sk to only the leading power in the quantum fields Xq,bq. To make the computation 
simple, we will carry this out separately for the saddle-point equation for A and bsp. 



A. Saddle point equation for A 



In order to compute 1/7V corrections to the saddle point equation for A (Eq. [29]), we write bd — > hsp + bd, and 
expand the action in Eq.[n]in powers of hd, bq and Ag. To achieve this we first integrate out the fermion fields in Eq.lS] 
to obtain 



Sk = -iNTrXn 



^mf ~ ^ bspT.c {bliTo + b*Tx) - {bclTo + bqT^) T,cbsp - {bclTo + bqT^) {b*ciTo + b*T^) 



'Xq idf — Ac 



N 



+2 b;) i^^g^ _a/ j \bq) + V'^J' ^^'^'p " '2\K{bd + bli) - 2Xdb.p {bq + b;) (57) 

Let us define, 

The solution to the above equation to leading order in A^ is 

Gmf = Gmf + XqGmfTxGmf = G,nf + XqSGmf 

where we define 

SGmf = GmfTxGmf 

Expanding Eq. [57]to quadratic order in the fluctuating fields bq^ci we get 
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{CmfiXq) (Ec (&*/To -I- b*qTx) + {bclTQ + hgT^) Ec 



GmfiXq){bdTo + bqrx)J:4b:,ro + blrx)] +2(6:, b^) U ^^t^^^A (bd 



+2Xq (^1 - - - 2Xqbsp {bd + b*i) - 2Xdbsp {bq + b;) 



(58) 
(59) 
(60) 



(61) 
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Collecting all terms upto quadratic order in the bosonic fields, we rewrite the action as below, 

Sk - -iNTr\nG;^){\q) + iNbspTr [G™/(A,) (S^ {KiTo + fo^r,) + (bciTo + b^r^) S^) 



-Xq idt - Ac/ 
idt - Ac; -A„ 



- n - -Xa( siia + sni^^ 



2(6ci bq 



N 



+2Xq ( 1 - — ) - 2Xqb% - 2Xqbsp {bd + b^i) - 2Ac/6.p (5, + bfj 



(62) 



n-' 



The above shows that the bosons due their interaction with fermions acquire the self-energies 11 = 

(5n(^'^^ = ( 2) ^) ) ' '^^'^ diagrams corresponding to 11, (jn'^'^^ are shown in Fig[2] (where the propagators 

are defined in Fig[T]). The bosonic self-energy 11 is, 



i'\ - (T''' [ToG.,nf{t,t')To^,{t' ,t)] Tr' [ToG™/(i,t')r,Sc(t',t)] 

^ ' ' 2 lTr'[r,G„,/(t,i')ToSc(t',i)] Tr' [T,G™/(t, t')r,Sc(t', t)] 



(63) 



where Tr' implies trace over only the Keldysh indices. Note that from causality the upper- left term in Eq. [63]is zero. 
Explicit expressions for 11 are given in Appendix [X] The other self-energies are. 



-iNbt„ 
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m«W(i,t') = ^^ 



(65) 



—iNh'^ 

<5n^W(i,0 = ^^ 



(66) 



which are of 0{1/N'^) and therefore will be dropped. The anomalous boson self-energies are the following 
— ?W/)2 



4 



(67) 



(68) 



c /ft 



(69) 



{{G^f^c)t^t' i^^f^^)t't + i^rnf^c)t,t' {^^f^f)t't + i^^f^c)t,t' i^^f^c)t,t + {G^f^c)t,t' (G^f^f) 



-lA ^A 



and are at least of 0{\/N). These will therefore not play a role in the C (;^) corrections to the saddle point equations, 

but will be important later, when we evaluate the conductance. The self-energies 5Il\ ' is also of 0{\/N) and wiU be 
dropped from further consideration. 
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FIG. 2: The bosonic self-energies corresponding to (a). 11; (b). 511^^^ and (c). ^n'^^ in text 



Other self-energies needed for computing corrections to the saddle point equations are 5Ilq and (5nJ. We find, 



_ -iiV /Tr' [ToSG„^fit,t')To^cit',t)] Tr' [ro(5G„/(t, i')TxSe(f, i)] 
2 \Tr' [TjGmfit, i')ToSc(i', *)] Tr' [r^SGrafit, t')T,^,{t' , t)] 

with SGmf defined in Eq. [501 Whereas, Su\^^ is (retaining terms upto 0{1/N)), 



(70) 
(71) 



jn«(t,t') = 



sn 



A{1) 



su 



R{1) 



(72) 



where 



« 2 



-iNbl 



G^f{t,t') (l]f<5G^^^E^)^, ^ + Gf„/t,t') (Sf <5G^„./Sf )^, ^ + Gl^(t,i') (S^^G^^Sf ),, J (73) 
G;:,/(t,t') (sf5G^/E^)^, ^ + Gif{t,t') (Ef JG^^Ef )^, ^ G^^(i,i') (E^<5Gj;,/Sf )^, 1 (74) 



We now integrate out the bosonic fields in the action Eq. [62l The ^^(A^) term cancels the last term because of the 
saddle point condition Eql^Il whereas the ©(Aj) term is also first order in the fluctuating bosonic fields 6g,c/- Thus 
on integrating out bq^ci, this term gives a term in the Keldysh action which is A^ and therefore does not affect the 
classical saddle point solutions. Following these steps we obtain. 



5 



K 



-2Xq ( 1 



iNTr\nG:;^^){Xq)+iTr\n 
N 



Dn^-n- A„ 



A„ ( m„ -t- 



- 2\qb' + 0{Xi) + Oil/N') 



(75) 



where i{babl) — ^I?o is the bare bosonic propagator. 

Now we may differentiate Eq[7S]with respect to A^ and set all quantum fields to zero in the resultant expression to 
obtain, 



2 1 



N 



+ 1 



NTr [G^j] 



iTr 



(76) 



Upto 0{1/N) only a subset of terms in Eq. [7S]need to be kept. Collecting these, 



2^1-1)- 2b% + ^NTr [G^,] - ^TrDl,. = f / ^ / 



dt f duj 
2^ 



(77) 



[D^^,4e)Gif{e + u;)G'^jie + Lo)+Dl^^,4e)G^jie + u;)G^^ie + uj) + D^^,,ie)Gif^^ {lo) 
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Substituting for the fermionic and bosonic Green's functions fEq. [24l [25l [33l IA121 IA141 IA15[) . we get 

A^f ^ r,/r 



i-bi 



a=L,R 



(78) 



r 



■K^Ti ^ r2 



dx- 



1 — iia- I'b 



-(Pa-Mb)/' 



1 ^ 1 



+ a; - ^i/r4 



ivrf ^ Tan 

-4 a,b=L,R 

We use the identity 



(1-^)^ 



1 + m^- 



r./r 



+ X- ^ii/TA 



1 



1 + mv - mx ^ ^— 



r 

r./r 



1 - ^iJTA 1 - 1^ 



1 + x - ^li/TA 



and for convenience introduce the foUowing short-hand, 



(79) 



Ta,b _ ^a^b 



D/Ta 



dx 



{pa-tJ.b)/TA 



2^ + "^E«T^ln|l+ J— 



(80) 



niT, ( X+{lLa-ll.i)lTA \ „^ , 



TjjT 



„n/T. f X+(pa-Hi)/TA \ 

J-i,.-M ""{x + mES^^l^ + T^ir 

Note that the functions /"'^ are infrared divergent. In equihbriuni (/i^ = /^i?, = 0) these have a logarithmic divergence, 
while out of equilibrium the 1°"'^ have a more severe l/x divergence. As shown in equilibrium by Read et al^ii^, in 
the computation of physical observables the 1°''' appear in such a way as to exactly cancel the divergences. For the 
out-of-equilibrium calculation as well we find an exact cancellation of divergences, so that all physical observables are 
well defined. 

In terms of the above symbols, Eg. [751 becomes 



— T^^-^ E \-mk'^'' + il + mv)P''' + Lr] (83) 



NT ^ Ta/T _ NTT 

a=L,R. -4 a,b=L,R 

The l.h.s of the above equation can be further arranged as follows by writing ep ^ Ta + P/N 



NT ^ rjr _NT ^ Ta/T p NT ^ Tg/T 

The 1/A^ correction to is carried out in the next subsection. Using the result for /3 derived there (Eq. llOip . the 
above equation is rewritten as 



AT Tg/T ^ AT r„/r _ / ^a,b \ NT_ Ta/T 

nep ^ l~ fia/ep ttTa ^ I ~ fia/TA [A^(l + ™y)^ ' j ttTa ^ {I - fia/TA^ 
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Substituting Eq. [85] into Eq. [83]we obtain the following expression for F upto 0{1/N) 



r 
f 

1 - 



mv 
1 + mv 



-y\- 

a,b \ 



-J^n — 



mv 



1 + mv 



(1 + mv)' 



r./r 



(1 - ^ijTAy 



It is convenient to introduce the following simplified notation 

a,b 
a.b 

r./r 



5'p=i,2,3 - XI n" 



(1 - fijTA)P 



Then, 



(86) 



1 m 



N l + mv 



(87) 
(88) 
(89) 

(90) 



r 1 + my 



m^ /mv m? jmv 



N \1 + mv ^ 1 + mv (1 + my)^ 
Note that in equilibrium when /i^ = /ij^ = 0, Eq. 1891 becomes 



p 



D/Ta ( 

dx ^ 



1 







X + mln (1 + x) 



1 

N 1 + mv' 



which is an expression that will appear later in the voltage expansion for physical observables. 

Thus the main result of this sub-section is Eq. [HI] which is the 1 /N correction to the level broadening. 



(91) 



(92) 



B. Saddle point equation for h 

In order to derive the 1 /N corrections to the saddle point Eq. [30l we set = in Eq. [6l write bd = bgp + hd , 6* 



bsp + b*i: bq = bq + bq,b* = bq + b*, and expand to quadratic order in the fluctuating fields bq^ci- Following this as 
before, we integrate out the fermions and the bosons and obtain an action Sxibq). To obtain the classical saddle 
point, we need -^j^-^li^^o- We will now follow the above steps. 
First we integrate out the electrons to obtain. 



Sk = -iNTrln 



Graf 



bspTo^cbqTx - bqT^T.cbspTo - 0(&^) - {bspT„ + bqT^) (6*;To + b*T^) 



(bclTQ + bqTx) Sc [bspTQ + bqTx) ~ {bclTQ + bqTx) Sc (^c/'^O + b*qTx)] 



-2 {b 



cl K) 



idt - Ac; 
idt - Xci 



) (^''') ~ ^■^'^'^'^p^i ~ 2Ac/6sp ipq + 6g) - 2Ac;6g (6c; + &c;) 



Expanding the above to quadratic order in the fluctuating fields we get 

Sk = im^pbqTv [G„/{r,Sc + SctJ] + iNb^bq {b^Tr [Gf„^Sf + G,^^I]^] + bq [G^^Ef + Gtj^^] ) 
+iNTr + bdbq) {G^fJ:^ + G^/E^) + {bqbsp + bqb*d) (G^;Ef + G;^/I]f )] 



'iNbq {bq 



'^{bli b*q) 



bspbq 



(93) 



'bd' 
. b„ , 



-2(&:; b*q 



2 (6c/ bq) 



-AXcibspbq - 2Xcibsp {bq + b*) - 2Ac/6g {bd + b*;) 



(94) 
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with n defined in Eq. [63l and the components of 511^ defined in Eq. [64l [65l [66l and those of SiP defined in Eq. 
and [Ml Moreover, the 611'^ are given by, 

An' i'\ - ziH ( ^'^ [(G™/(ScT, + T,Sc)G™/) (t, i')Sc(i', i)] Tr [(G™/(ScT, + t,Sc)G™/) (t, t)] , . 

' ^ " 2 ^^rr[T,(G,„/(S,T,+T,Sc)G™/)(t,t')Sc(t',i)] Tr [r, (G,„/(S,t, + r,Se)G„/) (t, t')T.Sc(i', *)] ^ ^ 

and ^Ilg = O ( and therefore will not play a role in the subsequent discussion. 
Integrating out the bosonic fields, and keeping terms upto 0{1/N ,bq) we get, 

Sk = iNbspbqTr [Gmf{Tx^c + ^cTx}] - ^Xclhspbq - ibspbqTr [DSU'q] 

Substituting for 511^, to 0{l/N), the above becomes, 

Sk = iNbspbqTr [Gm/jr^jEc + ^cTx}] - -iXcibspbq 
-Nbspbqj ^ I ^i?^(e)G,^,^(e + c.)E^(e + c.)Gf„^(e + o.)SfH 

Thus the saddle point equation for ep (obtained from = 0) reduces to 



(96) 



(97) 



eF-Eo + > In 



D 



(98) 



duj 



Far, 
b=L.,R r2 
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AT 



^ V J-^ In 1 1 ^ 



E 



duj- 



,„^ln|l ^ 



The above equation may be used to extract the 0{1/N) correction to the saddle point expression for the level 
energy ep- Writing 



Ta is given by Eq. [35l whereas from Eq. [98l we get 



^ ~ 1 + my r2 

a.h 



D/Ta 



dx 



I ^ 

\l-fia/TA l+X-llb/TAj 



iPa-tJ^b)/TA 



l + mE.^ln|l + ^ 



^J-^/TA I 



(99) 



(100) 



with my defined in Eq. 1371 
Using Eq. [5D] we may write 



m^^ m^ 
' ' " ^ 1 + my 



1 + my 



a, 6 



(101) 



Thus the two main results of this section is the 0{1/N) corrections to the level broadening (b'^p) and level position 
{Eq + A) which are given in Eq. [HI] and Eq. IIOII respectively. These results will be used in subsequent sections for the 
evaluation of various observables to O (■^). 



VI. EVALUATION OF np TO O 



In this section we will evaluate the local charge-density n^? to 0{1/N). np is given by 

N 



(102) 



Thus we need to evaluate to 0{1/N). For this we start by writing the Dyson equation for the fermionic Green's 
function correct to one loop, 



G/ — Gmf + Gmf^pGrnf 



(103) 



FIG. 3: Diagram contributing to the 1/N correction to np- 



where the second term in the above equation corresponds to the diagram in Fig[31 and the E^? are defined in Eq. lBli[B2l 
The Keldysh component of Eg. 11031 gives. 



We rewrite 
where 



NT 



Ta/r 



Y 



a—L,R 

—i ' 

-iN 



We use Eqns [OTl [Ml and [TUT] to correct f/r,ei. to 0{1/N) in Eq.dnHto obtain, 



my 



1 + my 



1 TO^ 

TV 1 + my 



Moreover using Eqn IBU IB21 one finds 



L°^' + (l + mv)r''' 



1 w? 

N l + mv 

a,b 

,3 



-E 



L 



a.h 



1 m mv \ ^ b 1 ™ 

TV 1 + my 



a,b 



N l + mv 



.a.b 



Adding Eqns [T06l [m] and [TT2] gives 

i + -y(- 

Comparing Eq. 11131 with Eq. |86]we may write the expressions for j: in the foUowing compact form, 



np — 



my 
1 + m,y 



Vmy my my 



■ mv 



1 + my 



(l + my)2 1 ^ (1 



r,/r 



/TaY 



r 1 m 

f = ^""^"iVTTmy^ 

a, 6 



E^' 



Upto C(l/A^), above may be rewritten as 



f = ^-nF-ii-np):!^yi-'' 



(104) 
(105) 

(106) 
(107) 

(108) 
(109) 

(110) 

(111) 
(112) 

(113) 



(114) 



(115) 
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The Ia,b contain the divergent terms. If A is an infrared cut-off, and defining V = |/iL ~ M-rL find, 

1 / r; 1 ,__^_ri 1 



r 1 

— ~ 1 — nj? 7T I ^ 7T 

r Nl + mv{l + mvf \ {1 - fiL/TA^ 



■In A 



1 



In -— 



(1 - ^iu/TAf 





r2 {^-flL/TA)il-^lR/TA) Ta (i - ^^/t^) (i - ^^^/t^) x 

The above expression will be useful in the next section when we study the bosonic correlation function 



(116) 



VII. BOSONIC CORRELATION FUNCTION: DECAY DUE TO CURRENT INDUCED 

DECOHERENCE 

The full bosonic correlation function (combining both saddle point and fluctuation corrections) is 



^^(t, t') = -^({6.p(^) + 5h{t), + 6b^it')}) = -2z- + i?^. (t, t') 
where -D^^. = ~i{{Sb{t) , 6b^ {t')}) , and its expression in frequency space is given in Eg. IA12I Using Eg. IA12I 



where 



(117) 



(118) 



7ab{t) = / dfle ' sgn {fl + fia ~ fJ-b) 



(119) 



where in the long time limit. 



7aa(i) 



~- In (AtTA] 



cos rtt 



(120) 



1-fj) (1 + m^)^ 



A is a cutoff introduced to take care of the infra-red divergences. This term, as we shall show will be canceled by the 
corresponding infrared divergence from Eq. 1911 
Similarly one finds (for V — — Hr\), 



ILR + IRL = 



Therefore the full bosonic correlation function is 



— cos nt + 2V — COS nt 



n 



V 



dfl 



02 



(121) 



-2i 



r 1 



r Nl + mv {l + mvf \ ^^{I-ixl/Ta 



■ In [MTa] - 



^ R 



1 



r2 (1 - ^jLlITa) (1 - ^iR|TA) Jv ^ 



°^dn 2rLrR 

—— cos i it - 



r2 (1 - ^,^/Ta) 
1 



■ In {MTa) (122) 



r2 (1 - f,L/TA) (1 - (ir/Ta) 



V I Y^cosm 
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(b) 



FIG. 4: Diagrams contributing to the susceptibility to 0{1/N). 



Combining the above with expression for F/F in Eg. 11161 one finds that the infrared divergences cancel to give, 



D^{i) = -2i{\-nF) 
2FLFfl 



1 - 



In {ITa) + 



F2 
^ R 



(1 - ^,j,/Ta 



■ In itTA) 



1 



F2 (1 - ^jLlITa) (1 - ^iR|TA) V Ta 



In— — h / cos lit 



V 



n 



1 



F2 (1 - (.l/Ta) (1 - fiRjTA) Ta Jo 



— (cosxtTA-l) 



(123) 



Note that the above expression is correct to 0{1/N). Therefore 1 — np needs to be computed only to the saddle 
point level for all terms except the first term in the square brackets. 
For long times Vt^^l, Eg. 11231 reduces to 



1 _ ££i intTA - — \ntTA - —Vt 
N N N 



where 



CLM. 



2 r2 
^ L,R 



1 



Cde 



[l + mvf r2 {l-^iL,RlTAf 

2TlTb. 1 



If one were to compute the correlation function to higher orders in , Eq. 11241 signals the following behavior 



(124) 

(125) 
(126) 



DK{t) ~ — 2« (1 — np) exp 



-l5:c,intr.-^i^. 

i=L,R 



(127) 



Thus the slow power-law decay in time in equilibrium of the bosonic correlator is replaced by a rapid exponential 
decay at non-zero voltages whose origin is current induced decoherence. Each of the exponents Ci^dec is consistent 
with what one might expect from nonequilibrium X-ray edge physics^^. 

The above decoherence rate appearing in the bosonic correlation function has consequences for physical observables 
such as the susceptibility and the conductance which we evaluate in subsequent sections. 



VIII. EVALUATION OF SUSCEPTIBILITY TO O (i) 



The spin response function is given by 



2 ^ 



duj 



27r 



(128) 
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ToO{l/N), from Eq.[T28l 



(129) 



where is the saddle point expression for the susceptibihty (diagram (a) in Fig [5]) with the level energy and 
broadening corrected to 0{1/N), while Xi'2 arise due to corrections to the electronic Green's functions to one-loop 
(diagram (6) in Fig[3|). In particular, 



m 

Xf ^^Y.{wBmfTr [{Gf„/SfG^^ + Sf G;^,/) + G^/S/G;^/}G^,^] 



(130) 
(131) 



In order to compute Xoi we use Eq. l46l and correct for f/F using Eq. [911 and correct for Ta using Eq. IIOII We 
find, 



Xo = Xsp 



1 



1 



N \ 1 + my " 1 + my 
Using Eq. IBll it is straightforward to show 

1 4 1 



m2 S3 



1 + mvSi (l + my)^ 



L1S2 — m I 



Xi 



R _ I ffVI 



J(J+ 1) 



Ta 3 iV V 1 + my 



■^3 = X, 



2 1 mig 
'"'3iV^ 



whereas 



X? = 2xf + 

9^^ilJ{J + l) 1 



TaFfc 

3 TV 1 + my ^ F2 

ab 



E 



(t^ 



(132) 



(133) 



(134) 



it^a.-IJ.b)/TA 



(x + m^,£fln|l + — £77^ 



r./F 



(1 + a; - fii/TA) 



Collecting all the terms together 



X^ = Xsp 



1 



, m^ 
k + 



N \ 1 + my 1 + my (1 + my) 



;LiS2 - 2 



^3 

S2I + my 



Li + 2mL3S2 - m^I + T 



where 



F 



^ F2 

a, 6 



■E 



D/Ta 



dx- 



{Pa-fJ.b)/TA 



r,;/F 



(1 + x - (i^/TaY 



Now it can be shown that 



(135) 



(136) 



where 



m 



1 + m,y 



-k-m^I + F 



1 + my 



-M 



Af = E A^""'' 

0,6 



Et 



F2 



D/Ta 



X+(Ala-Mb)/TA 

(l-Aia/TA)(l+an-A'b/TA) 

(^ + ^^,^r^ln|l + _£^|) 



1 + my 
5 



+ 2a: (1 - a^^/Ta) 



(1 + a; - ^Ji^/TAf{l - ^^^|TAf (1 + a: - ^Ji^/TA){l - ^i^/TA) 



(137) 



(138) 



(139) 



Therefore the static spin-susceptibility becomes, 



X^ = Xsp 



L2 



N \ 1 + my (1-1- my) 



;LiS2 - 2 



5,3 m2 
52 1 -I- my 



^3 

Li 2m— i - 



52 1 + my 



-Af 



(140) 
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with Xsp given in Eg. 1481 Eg. 11401 may be expanded in powers of pll,r/Ta- In particular in the Kondo hmit {jip 1 
or my ^ 1), Eg. 11401 is found to have the form, 



Xs 



np=i _ gV|J(J + 1) 



1 + 1.5 



^L^R f fJ-L - ^J■R 
p2 TO 



Ms 



vi=L,R 



^^^^ 

Tiny 



S2 



Ms 



ETO 

.i=L,R A 



Csi 



i=LM ^ 



{Cs3 ~ Csi) 



4.5 E^Efl / - fJ-R 
iV E2 I T° 



(Ul) 



The expressions for the Csi have been given in Appendix [Cl and may be evaluated numerically. The L^"? are defined 
in Eq. [Ml and 



jeq _ 
Jq — 



D/Ta 



dx 



(a; + mln(l + a;))^ {1 + xf 



(142) 



Let us assume that the chemical potential of the left lead fiL — V/2, and that for the right lead is fiR — —V/2. Let 
us define 



Cso = m {-Ll^ - Ll^ + 2Ll^ - m Jg^*) 



We now rewrite Eg. I141l as follows, 



(143) 



N 



1 / V 



N \ 2T\ 



Cs2 



1 fT,,-Tj 



N 



V 

2n 



p2 yj-O 



(C53 — Csi) 



N 

4.5 TlTr ( V 



N E2 \n 



N ) N 
2" 



\2Tl) 



Csi+ 
(144) 



where terms higher order than (r")^ have been dropped. Defining the Kondo temperature to 0{1/N) assi 



Tk = TI\1 



N 



(145) 



Eq. 11411 can be recast in the following universal form 



n^=i _ .gV|J(J + l) 
Xs — 



3Tk 



1 + 1.5 



.TlTr 



E2 \Tk 



V 



1 /T 

N 



V \ 



1 / V 



N \2Tk 



{Cs2 + — Csi) — — (4.5 + 3Cso + Cs3 — Csi) 



2Tk 
TlTr 



C 



si- 



V 



E2 VTa- 



(146) 



In the Kondo limit m ^ 1, the coefficients in the above equation take the following universal values Csi ^ 0.01, 
(Cs2 + Cs3 - Csi) ^ -0.005, (SCso + Cs3 - Csi) ^ -4.94. 



IX. EVALUATION OF THE SPECTRAL DENSITY AND CONDUCTANCE TO O {■^) 

The retarded Green's function, whose imaginary part gives the impurity spectral density is 
G%{t,t') = -^T{bUt)Ut)fUt')b-it'))-^{flit')^^^^^^ 
GL_{t, t')^^-' {t\ t) - GL+{t, t')^+-* (i', t) 



(147) 
(148) 



We are only interested in evaluating the imaginary part of Eq. 11481 which at leading order (saddle point level) is 
O {jf)- There are five diagrams that contribute to the above expression to O (775-) which are shown in Fig[5l Eor 
convenience we write 



Tn + Th + T, + 



(149) 
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FIG. 5: Diagrams needed for the computation of the impurity spectral density to 0{1/N'^) and hence the conductance to 



where Ti is the contribution from the z-th diagram and corresponds to 



Tain) 



1 



{Ta " nf (1 + mv)' 



1 + ^ - 



-L 



k 



N \ 1 + m,v 1 + mv (1 + my) 



1 + my _ 



duj 



r 



T,{n) = -Im 



Tdin) = 2-Im 



duj 



(150) 

(151) 
(152) 
(153) 



T,{n) = -Im 



doj 



+D^^,i^ + m^{~-u;)G^fi-co)]] (154) 



The above terms have been evaluated in Appendix [P] 

We now present results for the conductance for the case of symmetric couplings to the leads (r^ = Fj^) and 
Mi = V/2,HR = ~V/2. Defining 



Go — 



Ne'^ ATlTr / tt \ 2 / mo 



F2 



\n) \l + moJ 



(155) 
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and the functions 



DITa 



dx 



{x + m In (1 + x))^ 

2 



P = 
ti = 

<3 = 



(cc + m In (1 + x))^ 



E + 3) 



(i+-) 



(3x^+8z+6) 



x{x + 2) 



Q (a; + mln (1 + x)) 

1 f^/^- dx 

2 Jo {{x + m\n{l + x) f i'^ + x) 

-"/^■^ rfo; x{x + 2) X 

77 ^ i 7- '. 773 



1 
2 



1 



{i + xy 



{{x + m\n{l + x)f {l + xy (l + x) 
^/^-^ dx x{x + 2) 



10 {{x + m\n{l + xyf + 
we obtain the foUowing expression for the conductance, 



1 - 



l + x 



G{V) 
Go 
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as expected from the Friedel sum 



Note that in equilibrium, the above equation reduces to2^ G 
rule. 

The expression for the conductance in the Kondo regime may be obtained by setting too 3> 1 in the above equation. 
We find 
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Ne^ ATlTr /7r\2 
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with Tk defined in Eq. 11451 and the coefhcicnt Ggi given in Eq. IC7I In the Kondo limit, Ggi —2.77. 
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X. CONCLUSIONS 



In summary, we have presented results for the nonequilibrium infinite-U Anderson model using Keldysh functional 
integral methods. The approach has been to use 1/A^ as the small parameter in the theory which allows us to develop 
a systematic perturbation theory for the nonequilibrium problem. The results derived are valid for an applied voltage 
small as compared to the Kondo temperature when the effect of fluctuations are small. Physical quantities such as the 

impurity spin susceptibility and the conductance are calculated to O (^) ^ ' voltage expansions are found 

to show rich behavior by depending on different combinations of the couplings to the left and right leads such as: 

(^) 



— , '"pa" (^) • While terms of the first kind give rise to rectification type behavior, i.e. 

XsiV) 7^ Xs{—y) and G{V) ^ G(— F), the last term is associated with current induced decoherence as it arises due 



r2 
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to inelastic processes that can occur in an energy window V. This term is also found to cause the bosonic correlation 
function to decay rapidly in time. The approach developed in this paper is rather general, and therefore may be easily 
adaptable to a variety of out-of-equilibrium systems. 

An interesting question that arises is to what extent the results obtained in this paper are also valid for N — 2. It 
is known for equilibrium systems that a naive extrapolation of the results of large-N to iV = 2 when compared with 
exact Bethe-Ansatz results not only give incorrect numerical values of various quantities (such as the Wilson ratio 
and the zero-bias conductance), but also make qualitatively incorrect predictions for the temperature dependence of 
observables. Precisely how the extrapolation goes wrong has been discussed in Appendix [E] On the other hand, 
comparison with exact resultsii reveal that large-N works very well for > 4. However, one of the results of this 
paper has been the observation that G{—V) ^ G{V),xs{y) Xsi—V) for unequal couphng to the two leads. This 
asymmetry is rather generic and will exist whenever the system is away from particle-hole symmetry and therefore 
should be observed for the nonequilibrium N = 2 Anderson model away from the particle-hole symmetry point 
Eq — —-J- However, for a small voltage expansion of the conductance for N = 2 we do not expect the appearance 
of a linear in voltage term as found in Eq. [131 This is because the N — 2 case has a maximal conductance per 
channel of /h, and such a linear term would imply that the conductance can become larger than this value, which 

is unphysical. An asymmetry can very well appear at cubic order ( """^p^" (^) ^ small voltage expansion of 

G. Note that for A'' ^ 1, the conductance per channel is a small number of O (l/N"^). Therefore for this case a linear 
in voltage term in the conductance does not violate unitarity. 
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APPENDIX A: EVALUATION OF BOSONIC SELF-ENERGIES AND PROPAGATORS 



In this section we evaluate explicit expressions for H^^^^^ defined in Eq. [63l In particular 
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The above may be easily evaluated. We obtain. 
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Defining Xd + Eq — ep, and to 0{1/N), the above expressions simplify to 
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Similarly, to 0{l/N), H^^ is given by 
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The bosonic propagators may be evaluated from the Dyson equation 
To 0{l/N), (5n(i) does not contribute. So we have 



Evaluating the above, we get 

Re [Di^:,,m] 



Using the saddle point equation Eq. [35l the above becomes 
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Similarly, the imaginary part to 0(1/N) (where = Re[D^] + iIm[D^]) is 
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In the evaluation of the spectral-density, we also need the anomalous boson propagators Z?^ ^ 
0(i),wefind 

D^^, in) = in) = 2Dl,, in) sn^^'\n)D^^,, {-n) 

with ^-".^,^^(2) defined in Eqns. [m [BSland [Ml 



APPENDIX B: EVALUATION OF FERMIONIC SELF-ENERGY 



These are defined as 
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and are represented by the diagram in FigO 
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FIG. 6: Fermionic self-energy 

APPENDIX C: EXPANSION COEFFICIENTS IN THE EXPRESSION FOR THE SUSCEPTIBILITY 

AND THE CONDUCTANCE 

The expansion coefficients in Eq. 11411 are given by 

Cso = m (-^2' - Ll' + 2i^' - mJ^') (CI) 
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Whereas the coefficient Cgi appearing in the expression for the conductance in Eq. 11661 is 

Co. = -2toL^'' + ^ruLf - 6mi:^ - to^ f'^^ (^7) 

Jo (x + TOln(l + x)) U + 

APPENDIX D: COMPUTATION OF THE SPECTRAL-DENSITY 

In this section we give exphcit expressions for each of the five diagrams that contribute to the spectral density, 
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where Ti is the i-th diagram in Fig [5] We find 
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Note that the expression ID41ID5l and ID6l are logarithmically divergent. These divergences will be canceled by terms 
in diagrams (d) and (e), as we show below. In particular 
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Note that divergence in Eg. ID7I cancels the divergence in Eg. ID6I Moreover, we find 
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Note that Eq. IDIOI cancels the divergence in Eq. ID4I and ID5[ while the divergence in Eq. ID8I is canceled by that in 
Eq.[Di 



APPENDIX E: FAILURE OF EXTRAPOLATION TO N=2 FOR SYSTEMS IN EQUILIBRIUM 



Many N-fold degenerate magnetic impurity models besides being amenable to 1/N perturbative approaches, are 
also exactly solvable by Bethe-Ansatz. However, a comparison between results and exact solutions are not 
straightforward as the two approaches use different cut-off schemes (for a discussion on this see^i). Therefore the 
quantities that may be easily compared are universal quantities that are independent of the cut-off and Tk ■ One such 

quantitity is the Wilson ratio R — jj^^r^^r^ ^ , where xs is the impurity susceptibility and 7 = — §7^ is the specific 

heat coefficient. The exact solution of the Coqblin-Schrieffer Hamiltonian gives^^ 
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N - 1 
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Thus for N = 2, R = 2. On the other hand a perturbative expansion gives^^ R — 1 + Clearly, setting N — 2 
in this expression gives the rather incorrect result of i? = 1.5, showing that a naive extrapolation of the results of 
large N to the case of iV = 2 does not work. Another example is the value of the zero bias conductance through a 
TV-fold degenerate level in a quantum dot. The exact answer is 



G = N- 
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o TTHp 
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np being the average charge on the level which approaches the value ni? = 1 in the Kondo limit. Therefore when 

~ Tr and N — 2, the conductance in the Kondo limit reaches the maximum possible value of 2e^/h. For ^ 1, 

22 22 
a 1/A^ expansion gives G — j^j^, where again a naive substitution of A^ = 2 gives the incorrect result of G = 

The extrapolation besides giving incorrect numerical values, often does not capture the qualitative behavior of the 

temperature dependence of various observables. As an example let us consider the conductivity for the infinite-U 

Anderson model for a bulk system (rather than a quantum dot). If t~^(uj,T) = ilm [Gf^b] is the scattering rate due 

to the impurity, then the conductivity in a bulk geometry is cr''""'' ^ J du t{ll},T), whereas the conductance 

in a quantum-dot geometry (as has been considered in this paper) is G ^ J ^^f^^ T^^i^^T). The exact answer 
for a A^ = 2 bulk system is 
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where ct is a positive coefficient. On the other hand a 1/A" result for the temperature dependent conductivity for the 
infinite-U Anderson model ia^^ 



rbulk 
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Setting A^ = 2 in the above equation gives a qualitatively different result from Eq. IE3I as it predicts that the conduc- 
tivity will decrease with temperature (rather than increase). 

The above discussion shows that large-N results cannot be used to extrapolate to Af = 2. However, comparison 
with exact Bethe-Ansatz solutions-'^^ shows that large-N works well for A^ > 4. 
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